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Abstract 

The shape-dependent universality of the excess percolation cluster number and 
cross-configuration probability on a torus is discussed. Besides the aspect ratio of 
the torus, the universality class depends upon the twist in the periodic boundary 
conditions, which for example are generally introduced when triangular lattices are 
used in simulations. 



1 Introduction 



In standard random percolation, the universality of the critical exponents is 
quite familiar. The values of these universal quantities depend only upon the 
dimensionality of the system. This is also the case for the amplitude ratios of 
quantities such as the mean cluster size and correlation length that are defined 
essentially for infinite systems. Even though one might use finite systems to 
measure such quantities, in the limit of the system length scale L — > oo, the 
shape of those finite systems becomes irrelevant as L exceeds the correlation 
length 

However, if a system is simulated exactly at p c where £ is in principle infinite, 
or if one takes the scaling limit where p approaches p c as L is increased such 
that £ remains a fixed fraction of L, then the shape of the system (or the 
sequence of systems) becomes relevant. This leads to another and larger class 
of universal quantities, whose values depends upon the shape of the boundary 
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of the system and the boundary conditions, as well as the dimensionality. In 
spite of the resulting proliferation of classes, one still considers these quantities 
to be universal, because they remain independent of the "microscopic" features 
of the system (lattice type or continuum model used, whether site or bond 
percolation, etc.). 

Examples of shape-dependent universal properties go back to the work of 
Privman and Fisher [1] concerning Ising models on a torus. In percolation, the 
various crossing problems considered by Langlands et al. [2,3], Aizenman [4], 
Cardy [5], Ziff [6], Hovi and Aharony [7], Hu et al. [8,9], Watts [10] and Pinson 
[11] show this type of universality. These crossing quantites are all defined in 
terms of macroscopic features of the system, so universality is clear. Recently, 
Aharony and Stauffer [12] have examined the shape dependence of critical 
universal ratios such as L~ d S/ P^ where S is the mean size (in sites) of finite 
clusters and P^ is the probability a site belongs to the "infinite" cluster (one 
spanning the system, by some consistent definition). Although thess quantities 
are defined on a microscopic (site) level, the microscopic aspect cancels out by 
scaling theory and the ratio is universal — but again system-shape dependent. 
Similar quantities for the Ising model were studied by Kamieniarz and Blote 
[13] and Miiller [14]. 

Here we discuss two shape-dependent universal quantities: the excess clus- 
ter number and the cross-configuration probability. We also focus on twisted 
boundaries and the subtleties of the triangular lattice. 



2 The excess number of clusters 



Exactly at the critical threshold, the number of clusters per site or per unit 
area is a finite non-universal constant n c , whose values for site (S) and bond 
(B) percolation on square (SQ) and triangular (TR) lattices were examined 
in detail in [15]. Here we quote those results in terms of the number per unit 
area, taking the lattice bond length to be unity: nf -5 ^ = 0.0275981(3) and 
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which evaluates simply to [15] 
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Likewise, Baxter, Temperley, and Ashley's [17] integral expression for n B TR 
evaluates to (35/6 - 2/p^~ TR )y/3 = 0.129 146 645 . . . on a per-unit-area basis, 
where p B ~ TR = 2sin(7r/18). Here, clusters in bond percolation are character- 
ized by the number of wet sites they contain, with isolated sites corresponding 
to clusters of size one. 

While n c is non-universal, its finite-size correction is universal. For an L x V 
critical system with periodic boundary conditions on all sides, the total number 
of clusters N behaves as N(L, V) = n c LL' '+b+0(LL')~ 1 where 6 is a universal 
function of r = L'/L [15]. That is, 

b(r) = lim [N(L, L') - n c LL'} (3) 



with V = rL, and b represents the excess number of clusters over what one 
would expect from the bulk density. As such, it reflects the large clusters of the 
system, which is the basis of its universality. Aharony and Stauffer [12] showed 
that the universality of b follows directly from the arguments of Privman and 
Fisher [1] applied to percolation, implying that b is precisely the value of the 
free energy scaling function at p c and h = 0. Note that this "free energy" is 
J2s n s e ~ hs , different from the free energy given below. 

In [15], 6(1) was numerically found to equal 0.884 for both site and bond 
percolation on a SQ lattice, demonstrating universality. It was also found that 
6(2) = 0.991, 6(4) = 1.512, and for large r (systems of very high aspect ratio), 
6(r) ~ br with 6 = 0.3608. Here 6 is the excess number per unit length along 
an infinite periodic strip or cylinder. Periodic b. c. are essential in this problem 
to eliminate boundary effects which would otherwise overwhelm 6. 

In [18] it was shown that 6 can be found explicitly from exact coulomb-gas 
results. In the Fortuin-Kastelyn representation, the partition function of the 
Potts model at criticality is Z = Q Nc+Nb / 2 where Nc is the number of 
clusters and Nb the number of bonds, giving bond percolation in the limit 
of Q = 1. Here the free energy is F = InZ, and (N c + N B /2) follows from 
dF/dQ at Q — 1. Using the Potts model partition function (universal part) of 
Di Francesco et al. [19], [18] obtained 



h (r) = + g 5 / 4 ( 2v / 3r - \) + q 2 (V3r - 1) + g 5 / 48 + 2g 53 / 48 

_ g 23/16 + g 77/48 + (4) 

where q = e~ 2nr . This result yields 6(1) = 0.883 576 308 . . ., 6(2) = 0.991 781 515 .. . 
6(4) = 1.516 324 734..., and 6 = 5^/24, consistent with measurements of 
[15]. The result for 6 also follows directly from the work of Blote, Cardy and 
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Fig. 1. Two representations of a triangular lattice on a square array, yielding (a) a 
rectangular boundary of aspect ratio a/3/2 with no twist, and (b) the same rectan- 
gular boundary but with a twist of 1/2. 

Nightingale [20] on cylindrical systems. Fluctuations and higher-order cumu- 
lants are also discussed in [18], and b in 3d is discussed in [21]. 



3 Confirmation of universality using triangular lattices 

Demonstrating universality of b(r) by choosing only B-SQ and S-SQ as in [15] 
may not be completely convincing, and one would like to compare, say, a TR 
and SQ lattice. Periodic b. c. must be retained. Two obvious ways to represent 
the TR lattice on the square array of a computer program are shown in Fig 1. 
In (a), the TR lattice topology is created by choosing every other site of the SQ 
lattice. Taking an L x 2L boundary on the SQ lattice as shown, the effective 
boundary on the TR lattice becomes a rectangle with r = y/3/2. Periodic b. 
c. on the underlying SQ lattice results in normal (untwisted) periodic b. c. on 
the TR lattice also. The measured b = 0.887 for this system agrees completely 
with the theoretical b(y/3/2) = 0.887373 266 for a rectangular system with 
aspect ratio y/3/2. 

The second obvious way to represent the TR lattice — by far the most common 
one — is shown in Fig. 1(b). Diagonals are simply added to the SQ lattice, and 
the periodic b. c. are applied to the squared-off lattice as is. Making this into 
a proper TR lattice, the system becomes a 1 x 1 60° rhombus, and shifting 
around the triangle as shown in Fig. lb demonstrates that it is effectively 
a rectangular boundary with r = y/3/2, but with a "twist" t — \ in the 
periodic b. c, meaning that the x coordinates are shifted by a fraction t of the 
total length when wrapping around in the vertical direction. For this system, 
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Fig. 2. Excess cluster number b(r,t) vs. twist t for 14 x 16 simulation (triangles), 
theory for r = y/S/2 (solid line), 16 x 16 simulation (squares) and theory for r = 1 
(broken line). 

simulations gave b(y/Z/2, 1/2) = 0.878 [15], less than b(y/Z/2, 0) = 0.8874 and 
indeed less than the minimum untwisted rectangle, 6(1, 0) = 0.8836, where 
now we write 6 = 6(r, t) . 

Here we demonstrate the universality of b(\/3/2, 1/2) by studying a SQ lat- 
tice with r (necessarily rational) close to y/3/2 ~ 0.866, and comparing the 
results to the above TR lattice measurement. On the SQ lattice, we consid- 
ered a system of size 14 x 16, where r = 0.875; the measured values of 6 for 
t — 0, 1/8, ... , 1 are shown in Fig. 2. At twist |, the value of 6 is very close 
to the result 0.878 found on the TR lattice, [15] demonstrating the univer- 
sality between these two lattices. Note that, to find b(y/3/2, 1/2) on the SQ 
lattice to high precision, one would have to consider different size systems and 
extrapolate to oo, and different rational r to interpolate to r = y / 3/2. Results 
for a square system 16 x 16 are also shown in Fig. 2. 

We have generalized the theoretical methods described in [18] to find 6(r, t) 
from the partition function of [19]. The parameter r becomes t + ir, and 
the results, which are rather involved, yield the solid curves in Fig. 2. The 
discrepancy with the numerical values can be attributed to the small system 
size of these simulations. 



4 Symmetries on a torus with a twist 

The torus with a twist has various topological symmetries that apply to any 
shape-dependent universal quantity u(r,t), which includes b(r,t). We consider 
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a rectangular boundary with base 1 and height r, with a horizontal twist t 
in the periodic b. c. (Note that having twists in two directions leads to a 
non-uniform system, so we don't consider it.) u(r,t) satisfies the obvious 
symmetries of reflection 

u(r, t) = u(r, —t) (5) 

and periodicity in the t direction 

u(r,t)=u(r,l + t) (6) 

Another symmetry follows from the observation that the same rhombus can 
be made into a rectangle in two different ways, leading to: 

M(r ' t)=M (^'^) (7) 
Another construction shows that when t — 1/n where n is an integer, 

















\n z r 





(8) 



which also follows from Eqs. (5-7). On the complex r = t + ir plane, (7) corre- 
sponds to t — > l/r while (6) corresponds to r — > r + 1. These transformations 
generate the modular group, and functions invariant under them are called 
modular. Thus, b(r, t) must necessarily be a modular function. However, the 
explicit expression for b does not display that modularity clearly. 

Besides the excess number, another universal quantity on a torus is the cross- 
configuration probability 7r + (r, t), which can be expressed in a quite compact 
form. Using the results of [19], Pinson [11] has shown n + (r,t) = |[Z c (8/3) — 
Z c (2/3)], where 



Z c (h) = -J^- £ exp {- — K 2 + n 2 (r 2 + t 2 ) - 2tnn']\ , (9) 



rj(q) is the Dedekind eta function and q = e - 27r ( r ~ lt ) , it can be easily verified 
that this function satisfies the modular symmetries. For an untwisted torus 
(9) reduces to 

3 y(f)y(j:)-2y(fM£) 
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where fj(r) = i](e~ 2nr ) and tp(r) = i^e -7 ™") is the Jacobi theta function in 
Ramanujan's notation. The symmetry r — > 1/r is apparent. For all rational r, 
7r + (r, 0) is an algebraic number; for example, for r = 1, one can show 

Ml, 0) = 2 ^+^|- (' + ») + = Q.309526275 . . . (11) 

where a = 1 + \/3 and b = y/2 — 3 1 / 4 , using various results for theta functions. 
According to the symmetries above, this same value 7r+(l,0) = 0.3095... 
applies to (r,t) = (1/2,1/2) = (1/5,3/5) = (1/10,3/10) = (1/13,5/13) = 
(1/17, 13/17) = (1/26, 5/26) = (1/29, 12/29) = (1/34, 13/34) = (1/37, 6/37) = 
(1/41,9/41) = (1/50, 14/50) and infinitely many other systems, just as 6(1,0) = 
0.8836 . . . applies to all these systems. 

We have also measured n + in the 14 x 16 and 16 x 16 systems, with various 
t. Whether crossing occurs can be found using an indicator function, such as 
/ = Nc — Nc> + Nb — Ns, where Nc is the number of clusters, Nc is the 
number of dual-lattice clusters, Nb is the number of bonds, and N s is the 
number of sites. When 1 = 1, there is a cross-configuration on the lattice, 
when I = —1 there is a cross configuration on the dual-lattice (these two 
events are clearly mutually exclusive), and when 1 = there is neither. In 
the latter case, there will necessarily be at least one wrap-around cluster on 
the lattice or dual-lattice, or a spiral. Another indicator function can be made 
using the number of hulls in the system [11]. 

In Fig. 3 we show the measured n + (r,t) and comparison with predictions 
of Pinson's formula. The small deviations are presumably due to finite-size 
effects which should disappear when larger systems are measured and an 
extrapolation is made to infinity, as we have verified for t = 0. Note that 
7r + (v / 3/2, 1/2) = 0.316 053 413... is at a local and apparently global maxi- 
mum of Pinson's formula, although again by equations (5-7) an infinite num- 
ber of other points on the (r,t) plane have the identical value of n + , such as 
(V3/6,l/2) and (y/3/U, 5/14). 

A plot of Pinson's formula as a function of t for different r shows that for 
small r it begins to develop oscillations. This is because of the tendency to 
create spiral rather than cross configurations for small aspect ratios. 



5 The meaning of b 

In [15], it was suggested that b relates to the number of "spanning" clus- 
ters, since these are essentially the cause of the excess. However, Hu [22] has 
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Fig. 3. Pinson's number 7r+(r, t) vs. twist t; legend same as Fig. 2. 

shown that they are not numerically identical, using one particular definition 
of spanning. Here we elaborate on this point. 

We consider large systems with periodic b. c. and ask for the number of clusters 
Nr(£ m > £) whose maximum dimension £ m in the x or y direction exceeds some 
value £. Since s ~ £® where D is the fractal dimension and n s ~ s~ T , it follows 
that Nr(C > £) ~ £^ l ~ T ) D = t d or £~ 2 in 2d. We have measured this quantity 
for square Lx L and rectangular Lx2L systems for various sizes, and find that 
for an intermediate range in £ the expected universal £~ 2 behavior is followed: 



where A = L 2 (square) or 2L 2 (rectangular system) and C = 0.116. C is 
a universal measure of the size distribution, dependent only upon the rule 
of what constitutes the length scale £. (One could just as well use maximum 
diameter, radius of gyration, etc., each of which would lead to a different C. 
Thus C too is a "shape" -dependent universal quantity.) Eq. (12) implies that 
C is just the density in an infinite system of clusters of minimum dimension 
£, on the length scale of £. It is the universal analog of n c . 

The data for Nr(£ m > £) deviates from Eq. (12) at both the large and small 
size limits. For small £, the deviation is due to lattice-level effects; at the limit 
£ — 1, Nr(£ m > 1) is just n c A, which is clearly non-universal. For £ near the 
maximum, the deviation is due to the influence of the boundary and is related 
to the value of b. 

According to our definition, b = (the actual number of clusters) — (the ex- 
pected number of clusters using the bulk density). Therefore, using a lower 




(12) 
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length scale of £ that is in the scaling region, we have 



b = Nr(£ m >£)-CA/£ : 



a 



(13) 



Evidently, in using this formula, £ can be taken right up to the minimum 
dimension of the system, L. We find that the number in each size range, 



N t = Nr(£/2 < £ m < £) = Nr(C > - Nr(C > I) follows N e/2 /N e = 4 



within ±0.01 right up to £ = L. This implies that 6 can be found by applying 
(13) with I = L. Using the numerical data for Nr(£ m > L) for r = 1 and 
r = 2, we find from (13) 



compared with the actual vales 0.884 and 0.991 respectively. The small short- 
fall may be due to the need to take £ somewhat smaller than L in (13), or to 
statistical errors. We are investigating this point further. 

Thus, we have the meaning of b (taking I — L): the number of clusters in a 
system of area A whose extent is larger than the minimum system dimension 
L, minus the expected number predicted by the bulk density, CA/L 2 . 



6 Conclusions 

We have numerically demonstrated the universality of 6(v / 3/2,l/2) on both 
the SQ and TR lattices. We have shown that b is related to the average 
number of clusters of length scale greater or equal to L, if one subtracts off the 
contribution of the universal size distribution characterized by the universal 
constant C. 

It appears that Ng follows a universal behavior close to, and perhaps right up 
to £ = L. On a finite system with periodic b. c, the probability of growing a 
cluster of a certain size is identical to its probability on an infinite system, as 
long as the cluster is small enough that it doesn't touch itself after wrapping 
around the boundaries. This result, however, says that the number of clusters 
of a certain size range is also substantially unaffected by the system finiteness, 
even though the boundaries should, it seems, influence the statistics of such 
clusters when their combined size is large enough that they touch when wrap- 
ping around the boundary. Further work needs to be done to understand this 
behavior. 

The TR lattice constructed as in Fig. lb gives extremal values of b and 7r + and 



6(1,0) = 0.990 - 0.116 = 0.87 
6(2,0)= 1.214-2(0.116) = 0.98 



(14) 
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may in fact be the best periodic system to use for many percolation problems 
(as well as other types of lattice simulations). This is because, when the lxl 
60° rhombus is used to tile the plane, it leads to a triangular array of repeated 
patterns which has the most space between each repeated element of any 
regular array. 

While tori with "twists" add a rich extra degree of freedom, they are by no 
means the only systems for which b can be calculated. What is needed is a 
system that is effectively a closed surface. One can transform the rectangular 
basis of the torus to other shapes by conformal transformation, such as to an 
annulus, and then apply the transformed periodic b. c. to the problem. (Here, 
the curved boundaries suggest using a continuum form of percolation, as in 
[23]). A simple closed surface like the surface of a sphere can also be used for 
the system. Each of these systems will have its own characteristic value of b 
and other shape-dependent universal quantities. 
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